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Preliminary  Results  of  Low  Energy  Sputter  Yields  of  Boron 
Nitride  due  to  Xenon  Ion  Bombardment  (Preprint) 

James  L.  Topper*,  Binyamin  Rubin1,  Cody  C.  Farnell^,  Azer  P.  Yalin§ 

Colorado  State  University,  Fort  Collins,  CO,  80523 

We  describe  the  development  of  an  experimental  apparatus  geared  towards 
measurement  of  boron  nitride  sputtering  by  low  energy  ions.  A  four-grid  system  is  used  to 
achieve  a  collimated  beam  at  low  energy  (<100  eV).  A  weight  loss  approach  is  used  to 
measure  total  sputter  yields  and  a  quartz  crystal  microbalance  (QCM)  is  used  to  measure 
differential  sputter  yield  profiles  of  condensable  components.  Integration  of  the  QCM 
profiles  also  gives  total  sputter  yields  of  condensable  components.  We  report  initial  results  of 
total  and  differential  sputter  yield  measurements  of  three  grades  of  boron  nitride  due  to 
bombardment  by  xenon  ions  with  ion  energies  in  the  range  of  60-250  eV  and  at  ion  incidence 
angles  of  0, 15,  and  45  degrees  from  normal.  Comparison  with  past  measurement  results  are 
made  where  possible. 


I.  Nomenclature 


A, 

= 

sensor  area  of  QCM 

E 

= 

beam  ion  energy 

E* 

= 

characteristic  energy  to  describe  sputtering  profile 

J B,  avg 

= 

time-averaged  current  of  bombarding  particles  (ions  and  energetic  neutrals) 

Mt 

= 

molar  mass  of  species  i 

R 

= 

mass  accumulation  rate 

Y 

= 

total  sputter  yield 

a 

= 

polar  angle  of  sputtered  particle  (from  target  normal) 

P 

= 

ion  incidence  angle  relative  to  target  normal 

Am 

= 

mass  loss  due  to  ion  beam  exposure 

0 

= 

azimuthal  angle  of  sputtered  particle 

P 

= 

density  of  target  material 

f qcm 

= 

distance  from  the  target  center  to  the  QCM 

Ymz 

= 

modified  Zhang  differential  sputter  yield 

y 

= 

volumetric  differential  sputter  yield 

II.  Introduction 

Ion  sputtering  is  the  process  in  which  atoms  (and  molecules,  clusters,  or  ions)  are  ejected  from  the  surface  of  a 
material  due  to  bombarding  incident  ions1.  Details  of  the  sputtering  process  are  of  interest  from  both  theoretical 
and  applied  viewpoints.  Our  primary  interest  is  to  better  understand  the  role  of  sputtering  in  electric  propulsion  (EP) 
thrusters  used  for  satellite  and  space  exploration2'10.  More  specifically,  we  are  interested  in  the  sputtering  of  boron 
nitride  (BN)  because  of  its  widespread  use  as  an  insulator  material  in  the  acceleration  channel  of  stationary  plasma 
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thrusters  (SPTs).  Erosion  of  the  insulator  channel  is  an  important  thruster  life-limiting  mechanism,  and  deposition 
of  the  sputtered  BN  can  contaminate  spacecraft  surfaces  (e.g.  solar  panels  or  thermal  control  surfaces). 

Owing  to  relatively  long  lifetimes  (5-10+  years)  of  EP  thruster  devices  and  the  complexity  and  expense  of 
experimental  tests,  effects  of  sputter  erosion  and  deposition  are  generally  studied  with  numeric  codes.  For  erosion 
studies  (lifetime),  one  aims  to  compute  the  amount  of  surface  erosion  due  to  the  bombarding  ions.  The  modeling 
requires  knowledge  of  the  total  sputter  yields  (Y)  of  the  eroding  materials  of  interest  at  the  ion  conditions  of  interest. 
Modeling  of  deposition  additionally  requires  knowledge  of  the  differential  (angular)  sputter  yields  ( y(a,<?> ))  in  order 
to  track  the  trajectories  of  sputtered  particles.  Total  and  differential  sputter  yield  profiles  have  been  measured  with  a 
multitude  of  techniques,  a  partial  list  of  which  includes  weight  loss11,  collector  plates12'13,  mass  spectrometry14, 
quartz  crystal  microbalance  '  ’  "  ,  Rutherford  backscattering  "  ,  radioactive  tracers  ,  and  cavity  ring-down 
spectroscopy20. 

Despite  the  importance  of  BN  erosion  there  is  a  lack  of  fundamental  sputtering  data  on  BN.  In  particular, 
numerical  modeling  of  thruster  erosion21'23  shows  that  the  ions  most  critical  to  the  erosion  process  have  relatively 
low  energy  (<-100  eV).  Given  the  low  associated  sputter  yields,  such  measurements  strain  detection  limits  of  many 
measurement  systems,  making  them  especially  challenging.  The  goal  of  the  present  submission  is  to  contribute 
towards  filling  this  gap.  In  this  work  we  detail  development  of  an  experimental  configuration  for  low  energy  BN 
sputter  measurements.  We  report  initial  results  for  several  grades  of  BN  at  ion  energies  down  to  and  below  100  eV. 
The  results  presented  here  represent  the  start  point  of  a  comprehensive  study  that  is  currently  in  progress  and 
encompasses  a  broader  range  of  sputtering  conditions  and  includes  temperature  variation  studies.  The  measurement 
approach  is  to  use  a  combination  of  weight-loss  and  QCM  deposition  sensor  and  builds  upon  our  previous  work5'7, 9' 
10.  In  Section  III  we  summarize  our  experimental  apparatus  including  sputtering  chamber,  ion  source,  and  weight 
loss  and  QCM  measurement  systems.  In  Section  IV  we  present  sputtering  results  for  bombardment  of  xenon  ions  on 
HBC,  HBR,  and  HP  grades  of  boron  nitride  including  comparison  with  other  measurements  (where  possible). 
Finally,  conclusions  are  given  in  Section  V. 


III.  Experimental 

A.  Overview  of  Sputter  Measurement  System 

Sputter  measurements  are  performed  using  a  quartz  crystal  microbalance  (QCM)  deposition  monitor  as  well  as  a 
weight  loss  approach.  An  improvement  relative  to  our  past  work  is  that  both  measurements  are  performed 
concurrently  in  the  same  facility.  The  experimental  apparatus  is  shown  in  Figure  1 .  The  system  has  been  previously 
described5'7,  910  so  that  in  this  subsection  we  give  an  overview  of  its  essential  features,  while  the  following 
subsections  detail  specific  aspects  and  recent  modifications.  The  ion  source  and  QCM  are  housed  within  a  0.125  m3 
stainless  steel  vacuum  chamber  (43  cm  ID  x  76  cm  long  main  section),  equipped  with  a  1500  liter/s  CTI-8  cryogenic 
pump.  The  chamber  base  pressure  was  5xl0'7  Torr  giving  a  working  pressure  of  approximately  0.6  to  lxlO'4  Torr. 
The  gridded  ion  source  has  been  specially  designed  for  low  energy  operation  and  is  described  below  in  more  detail. 

A  rotatable  target-mount  is  positioned  23  cm  downstream  of  the  ion  source.  Details  of  the  BN  targets  used  are 
presented  in  Section  IIIC.  We  have  incorporated  a  heating  plate  into  the  sample  holder  as  a  means  to  perform 
sample  bake-outs  and  to  control  sample  temperature  for  upcoming  temperature  dependence  studies.  The  new 
sample  holder  is  formed  from  a  6  x  6  x  1  inch  plate  of  stainless  steel  with  four  holes  drilled  into  a  6  x  1  inch  side. 
These  holes  are  evenly  spaced  and  each  penetrates  4  inches  into  the  plate  with  a  diameter  of  0.5  inches.  Into  each  of 
these  holes  we  place  a  cylindrical  cartridge  heater  capable  of  outputting  as  much  as  200  W  of  power  (per  cartridge). 
At  full  power,  the  plate  can  raise  the  sample  temperature  as  high  as  700  °C,  although  other  constraints  of  the  system 
will  limit  the  temperature  to  about  500  °C. 
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Figure  1.  Schematic  diagram  of  experimental  set-up 


B.  Four-Grid  Ion  Source 

We  have  developed  a  four-grid  ion  source  to  enable  collimated  beams  at  the  low  ion  energies  of  interest  to  this 
work.  Detail  on  the  design  and  manufacture  of  the  source  will  be  reported  in  a  future  submission.  Here,  we 
summarize  the  operating  conditions  and  ion  beam  characteristics.  The  source  is  designed  to  operate  on  xenon  gas 
(typical  flow  rate  of  0.5  seem)  and  is  based  on  a  plasma  discharge  chamber  operating  with  a  dual  filament  cathode. 

The  four-grid  ion  source  consists  of  a  screen  grid,  an  acceleration  grid,  a  focus  grid,  and  a  deceleration  grid.  The 
screen  grid  controls  the  beam  energy  and  is  maintained  positive  to  ground,  while  the  acceleration  grid  draws  the  ion 
beam  out  of  the  source  and  is  maintained  at  approximately  -600  V  relative  to  ground.  The  focus  grid  collimates  the 
beam,  and  is  held  negative  to  ground  generally  at  -125  V.  Finally,  the  deceleration  grid  helps  prevent  ion 
backstreaming,  and  is  generally  maintained  slightly  negative  (~2  V)  relative  to  ground.  The  beam  then  passes 
through  a  ground  screen  and  over  a  neutralizer  filament  kept  at  an  emission  current  of  roughly  twice  the  beam 
current.  The  filament  passes  directly  through  the  center  of  the  beam  to  maximize  its  efficiency. 

Beam  profiling  tests  were  performed  using  a  Faraday  probe  mounted  on  a  multi-axis  rotatable  base.  Profiles 
were  integrated  to  find  total  beam  current.  For  typical  operating  conditions  the  radius  encompassing  90%  of  the 
total  beam  current  was  - 2  cm.  Table  1  shows  extracted  total  beam  currents  at  tested  beam  energies.  We  were  able  to 
extract  high-current  and  well-collimated  beams  even  at  energies  as  low  as  30  eV. 


Table  1.  Beam  energies  and  currents  from  four-grid  ion  source. 


Beam  Energy  (eV) 

30 

40 

60 

80 

100 

150 

200 

250 

Beam  Current  (mA) 

0.97 

1.14 

1.46 

1.96 

1.84 

2.11 

2.94 

3.45 

During  sputter  measurements  we  record  (from  the  power  supplies)  the  ion  current  leaving  the  source. 
Determination  of  the  sputter  yield  requires  the  current  of  energetic  particles  (ions  and  fast  neutrals)  incident  on  the 
sputtering  target.  As  in  our  past  work9,11,  we  make  corrections  for  charge-exchange  and  scattering.  The  charge- 
exchange  beam  generates  fast  neutrals  which,  depending  on  scattering  angles,  may  bombard  the  target.  The 
resulting  correction  is  to  multiply  the  source  current  by  0.95  ±  0.05. 

C.  Boron  Nitride  Targets:  Surface  Charging  and  Moisture  Effects 

Test  results  reported  herein  are  for  HBC,  HBR,  and  HP  grades  of  Boron  Nitride  (BN).  Each  of  these  materials  is 
formed  by  hot-pressing  and  corresponds  to  the  graphite-like  allotrope  of  BN.  In  the  base  plane,  atoms  are  held 
together  by  strongly  directed  hexagonal  arrays  of  covalent  bonds,  resulting  in  impressive  electrical,  thermal,  and 
mechanical  properties.  The  HBC  and  HBR  materials  were  obtained  from  General  Electric's  (GE’s)  Advanced 
Ceramics.  Calcium  borate  is  used  as  binder  in  HBR,  however  no  binder  is  used  in  HBC.  The  two  grades  have 
generally  similar  properties  though  with  some  differences.  For  example,  HBR  has  higher  thermal  expansion,  higher 


3 

American  Institute  of  Aeronautics  and  Astronautics 


moisture  absorption,  and  higher  volume  resistivity  at  elevated  temperatures.  More  detail  on  the  materials  can  be 
obtained  from  GE  datasheets.  The  HP  grade  is  obtained  from  Saint-Gobain  and  uses  calcium  borate  as  a  binder. 

Our  past  measurements  of  insulator  materials,  including  boron  nitride,  have  shown  effects  of  surface  charging 
and  the  importance  of  appropriate  neutralization.  Similar  surface  charging  effects  have  been  observed  by  Zhang  et 
al.24  and  Nikiporetz  et  al.25.  In  order  to  neutralize  the  surface  charge,  a  plasma  bridge  neutralizer  (PBN)  was  placed 
in  the  chamber  close  to  the  sample  being  sputtered.  Details  on  the  neutralization  scheme  have  been  presented  in  our 
past  work10.  Operating  conditions  of  the  PBN  include  an  emission  current  of  10-15  mA  and  a  Xe  mass  flow  rate  of 
0.5  seem.  The  PBN  was  biased  negatively  relative  to  ground  potential,  at  roughly  -15  V. 

Our  past  measurements  have  shown  that  HBR  BN  exhibits  effects  of  moisture  absorption  and  associated  mass 
change  at  relative  humidity  levels  of  less  than  or  equal  to  40%  at  room  temperature,  while  HBC  does  not10.  These 
effects  are  expected  based  on  moisture  absorption  information  provided  in  the  material  datasheets.  Similar  effects 
were  reported  by  Gamier26.  HBR  BN  can  maintain  appreciable  amounts  of  moisture  even  after  storage  in  a  dry 
environment.  The  exact  mechanism  of  absorption  is  unclear  but  it  is  postulated  that  the  surface  undergoes  chemical 
or  compositional  changes  related  to  the  calcium  borate  binder.  Empirically,  we  find  that  (for  the  resolution  of  our 
measurements)  the  mass  buildup  levels  off  after  approximately  one  hour  after  venting  the  chamber  and  exposing  the 
targets  to  the  atmosphere.  To  accommodate  this,  the  approach  we  follow  is  to  perform  the  before-  and  after¬ 
measurements  in  the  heavy  (moist)  state  by  waiting  one  hour  after  removal  prior  to  measuring  the  mass.  The 
method  has  been  experimentally  validated.  In  our  previous  work,  we  have  discussed  baking  the  samples  for  ~0.5 
hours  to  remove  absorbed  moisture.  In  the  present  study,  the  pre-test  bake-out  was  eliminated.  Instead,  the  samples 
were  heated  under  the  chamber  filaments  prior  to  sputtering  and  only  allowed  to  reach  a  thermal  steady- state 
condition  (-20  minutes).  Since  the  mechanism  behind  the  moisture  absorption  is  not  completely  understood,  it  was 
decided  that  baking  the  samples  at  higher  temperature  could  potentially  cause  changes  in  both  surface  structure  of 
the  sample  and  the  chemical  composition  of  the  calcium  borate  binder,  which  could  alter  sputter  properties. 

Similar  issues  arise  for  HP  BN  and  limit  our  ability  to  do  sputter  measurements  by  the  weight  loss  method.  The 
HP  BN  also  contains  the  calcium  borate  binder  (and  in  higher  mass  percentage)  so  that  moisture  effects  may  be 
expected.  In  comparison  to  HBR  BN,  we  found  the  HP  mass  to  increase  at  a  higher  rate  and  for  a  longer  duration 
(more  than  24  hours).  As  means  of  illustration,  Fig.  2  shows  mass  buildup  on  BN  samples  after  removing  from  the 
vacuum  chamber.  The  curves  show  mass  buildup  beginning  after  90  minutes;  mass  build-up  in  the  first  90  minutes 
is  at  an  even  higher  rate.  In  summary,  it  was  found  that  the  large  mass  buildup  (and  variation  of  individual  samples) 
relative  to  the  total  before/after  mass  change  of  a  typical  test  (-2-10  mg)  precludes  weight  loss  measurements  of  HP 
BN  in  our  current  setup.  See  also  discussion  of  humidity  effects  in  Section  IVA.  Such  effects  do  not  influence  QCM 
measurements  so  total  sputter  yields  for  HP  BN  can  be  inferred  from  QCM  results. 


Time  (minutes) 


Figure  2.  Mass  accumulation  as  a  function  of  time  from  six  different  HP  BN  weighing  sessions.  Time  zero 
corresponds  to  90  minutes  after  removal  from  vacuum  chamber. 
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D.  Definition  of  Angles 


The  angles  used  to  describe  the  direction  of  ion  incidence  and  the  ejections  angles  of  sputtered  particles  are 
shown  in  Fig.  1.  We  define  as  follows:  /?  is  the  incidence  angle  of  bombarding  ions  measured  relative  to  the  surface 
normal  (/?=0  for  normal  incidence),  a  is  the  ejection  polar  angle  of  sputtered  atoms  measured  relative  to  the  surface 
normal,  and  ^is  the  ejection  azimuthal  angle  of  the  sputtered  atoms  measured  in  the  plane  of  the  target  surface 
(defined  so  that  ^=0  is  in  the  forward  sputter  direction  i.e.  in  the  forward  direction  of  the  plane  containing  the 
surface  normal  and  the  incident  ion  directions). 


E.  QCM  Sensor  and  Signal  Analysis 

In  deposition  mode,  the  QCM  allows  determination  of  differential  sputter  yields  through  measurement  of  mass 
accumulation  (of  sputtered  particles)  on  its  surface.  For  condensable  components,  sticking  coefficients  are  assumed 
to  be  unity.  Note  that  sticking  coefficients  for  “new  layers”  and  very  thin  layers  (Angstrom  thicknesses)  may  be  less 
than  unity,  but  once  a  sufficient  layer  thickness  of  a  given  material  has  accumulated,  sticking  coefficients  for 
condensables  are  generally  unity27.  For  the  BN  materials,  the  sputtered  particles  may  consist  of  a  mix  of 
condensable  and  noncondensable  components;  this  issue  is  further  addressed  in  Section  IV. 

We  use  a  Sigma  Instrument  SQC-339  Deposition  Controller  that  reads  the  crystal  frequency  to  0.001  Hz  and  an 
RC-cut  quartz  crystal  as  opposed  to  the  more  conventional  AC-cut  crystal.  The  RC-cut  crystal  manufactured  by 
Tangidyne  Corporation  is  extremely  accurate  for  deposition  of  very  thin  films.  Increased  sensitivity  is  achieved  by 
adjusting  the  stress  coefficients  of  the  quartz  plate  using  advanced  fabrication  methods.  For  most  materials  of 
interest,  the  RC-cut  crystal  enables  measurement  of  differential  sputter  yields  at  ion  energies  below  100  eV. 

Since  quartz  crystal  resonance  frequency  is  extremely  sensitive  to  temperature  variation,  it  should  be  maintained 
at  constant  temperature  during  the  sputter  yield  measurements.  A  Poly  Science  9002  Programmable  Digital 
Temperature  Controller  with  circulating  bath  is  used  to  control  the  temperature  of  the  QCM.  The  water  is  circulating 
through  the  stainless  steel  body  of  the  QCM  housing,  while  the  temperature  of  the  water  in  the  bath  is  controlled  to 
0.01  °C.  As  the  crystal  is  moved  to  different  positions  during  the  measurement,  the  heat  flux  to  it  varies  due  to  the 
change  in  relative  position  of  the  QCM  and  the  heat  sources  in  the  system,  that  is,  the  ion  source  and  PBN. 
Therefore  although  the  temperature  of  the  water  stays  constant,  the  actual  crystal  temperature  is  different  at  different 
locations  (wangles).  A  K-type  thermocouple  wire  embedded  in  a  copper  holder  silver-soldered  to  the  back  of  the 
QCM  crystal  holder  is  used  to  monitor  QCM  temperature,  and  a  Lab  View  program  that  performs  data  logging 
during  differential  sputter  yield  measurements  has  been  modified  to  ensure  that  the  temperature  of  the  QCM  is  stable 
during  every  measurement.  When  the  QCM  is  moved  to  the  next  measurement  position,  the  program  monitors  its 
temperature  and  starts  the  sputter  yield  measurements  only  after  the  temperature  of  the  crystal  has  stabilized  to 
within  0.02  °C.  The  improved  temperature  stability  aids  in  high  signal-to-noise  ratio  sputter  measurements  as  are 
needed  for  low  ion  energies.  The  improvement  over  the  previous  system  is  largely  due  to  the  direct  monitoring  of 
the  QCM  temperature  (as  opposed  to  monitoring  of  the  temperature  of  water  in  the  circulating  bath).  The  criterion 
used  for  temperature  stability  control  can  be  adjusted  if  necessary  to  achieve  even  higher  thermal  stability,  however 
each  measurement  would  take  substantially  longer  to  acquire. 

For  a  given  incidence  angle  (obtained  by  tilting  the  target),  the  differential  sputtering  profile  is  obtained  by 
measuring  the  sputter  yield  over  two  chords  above  the  target:  </)=  0°/180°  and  (/)=  60°/240°  (where  the  latter,  by 
symmetry,  also  corresponds  to  (j)  =  120°/300°  for  azimuthally  symmetric  sputtering  profiles).  A  total  of  ~34 
positions  above  the  target  are  typically  sampled.  At  a  given  measurement  point  the  volumetric  differential  sputter 
yield,  y((X,(/)),  in  units  of  mm3/C/sr,  is  determined  using  equation  (1).  Note  that  because  we  study  multi-component 
materials  (for  which  the  sputtered  particles  may  comprise  various  atoms  or  molecules),  we  use  volumetric  units,  e.g. 
mm3/C/sr,  as  opposed  to  atomic  based  yields,  e.g.  atoms/ion.  In  Eqn.  (1),  R(a,</))  is  the  measured  mass  accumulation 
rate  (found  from  a  deposition  monitor  device),  p  is  the  density  of  target  material,  JB,avg  (C/s)  is  the  time-averaged 
current  of  bombarding  particles  (ions  and  energetic  neutrals)  incident  on  the  target,  rqcm  is  the  distance  from  the 
target  center  to  the  QCM  (17.4  cm),  and  As  is  the  QCM  sensor  area  (0.535  cm2).  The  quantity  AJrqcm 2  corresponds 
to  the  solid  angle  that  the  QCM  sensor  subtends  while  R(a,  (/))/pJB,avg  corresponds  to  the  volume  of  sputtered  material 
per  bombarding  charge.  Again,  it  is  important  to  emphasize  that  the  directly  measured  quantity  is  the  mass  buildup 
of  condensable  particles  on  the  QCM  and  the  volumetric  differential  sputter  yield  should  be  considered  in  this  way. 
(In  fact  the  “volume”  may  not  really  correspond  to  any  physically  observed  volume  since  it  corresponds  to  the 


5 

American  Institute  of  Aeronautics  and  Astronautics 


equivalent  volume  due  to  the  mass  of  the  deposited  condensable  material  if  one  uses  the  density  of  the  full  target 
material ;  of  course  conversion  to  the  deposited  mass  yield  simply  requires  removing  the  density  from  Eqn.  (1)) 

y{a,(f))  =  [/?(#,  (p)rqcm  B,avgA ]  (D 

In  principle,  because  of  the  finite  size  of  the  QCM  crystal  and  beam  spot  on  the  target,  each  of  our  measurements 
(i.e.  QCM  positions)  corresponds  to  a  (small)  range  of  polar  and  azimuthal  angles  joining  the  target  and  QCM.  We 
have  performed  a  simple  simulation  to  show  that  for  our  geometry  and  sputter  conditions  these  effects  are  negligible 
(worst  case  of  5%  error),  so  that  we  can  treat  our  target  and  QCM  as  points  (not  areas)  joined  by  a  single  vector5. 
We  also  note  that  the  viewing  angle  of  the  QCM  is  approximately  5.4  steradians  (cone  apex  angle  of  165°)  which 
corresponds  to  an  area  on  the  target  that  is  larger  than  the  irradiated  area  (also  a  necessary  requirement  for  equation 
(1)).  Experimental  uncertainty  and  error  bars  on  measured  differential  sputter  yields  are  found  by  estimating 
individual  contributions  and  combining  them.  Uncertainties  are  as  follows:  ±  15%  on  beam  current,  +  6%  on  radial 
distance  from  QCM  to  target,  and  +  0.1%  on  measurement  of  the  a  and  (j)  angles.  The  resultant  uncertainty  varies 
depending  on  measurement  conditions  but  typical  uncertainties  are  approximately  ±  16%.  For  simplicity  in  analysis, 
we  assume  that  the  uncertainty  in  total  yield  Y  and  E*  from  the  QCM  have  the  same  fractional  value. 

Analysis  and  fitting  of  differential  sputter  yield  profiles  requires  appropriate  functional  forms.  At  our  conditions, 
stopping  is  predominantly  due  to  elastic  (nuclear)  collisions  and  is  generally  in  the  linear  cascade  regime  (emitted 
particles  are  secondary  or  higher  generation  recoils)  or  single  knock-on  regime  (emitted  particles  are  primary 
recoils)1.  A  classical  theory  for  the  linear  cascade  regime  was  originally  developed  by  Sigmund28.  Independent  of 
ion  incidence  angle,  the  original  Sigmund  theory  predicts  sputtering  profiles  that  are  azimuthally  symmetric  and 
approximately  diffuse  in  shape,  corresponding  to  cosine-like  profiles  of  the  form  yoccos(dr)n  (n=  1  for  a  diffuse 
profile).  More  recent  experimental  and  numerical  studies  show  a  range  of  profile  shapes.  For  normally  incident 
ions  on  polycrystalline  and  amorphous  targets,  cosine-like  profiles  are  generally  observed  with  increasingly  under¬ 
cosine  shapes  as  ion  energy  is  lowered  and  increasingly  over-cosine  shapes  for  higher  ion  energies5'7,9'10,14,29'31.  For 
obliquely  incident  ions  at  relatively  high  ion  energy,  observed  profiles  also  tend  to  be  azimuthally  symmetric. 
However,  for  lower  ion  energies  the  measured  profiles  tend  to  be  asymmetric  with  increased  sputtering  in  the 
forward  direction5'7,9'10,15,19,29.  Similar  profiles  have  been  modeled  on  a  theoretical  basis32'34. 

As  a  means  to  describe  the  measured  differential  sputter  yield  profiles  we  use  expressions  from  Zhang35,  based 
on  work  from  Yamamura  33'34,  to  which  we  introduce  two  fit  parameters.  We  term  the  resulting  expressions  as 
Modified  Zhang  (MZ)6,14. 


y  mz 


cos  (or) 


l-,r  eos(/>’) 
E 


71 


1  E*  (  3 

1 — .  —  cos(/3)y(a)  +  —7rsin(j3)sin(a)cos( 
4  V  E  l  2 


(2a) 


^  _  3 sin  ( af  - 1  +  cos  ( < osf  (3 : sin  ( < a)2  + 1)  ^  j'  i  +  sin  ( a) 


sin  (dr) 


2 sin  (dr) 


1-  sin  (dr) 


(2b) 


where  yMZ  is  the  differential  sputter  yield,  Y  is  the  total  sputter  yield,  E  is  the  ion  energy,  F*  is  a  characteristic 
energy  describing  the  profile  shape,  and  the  angles  are  as  defined  above.  The  approach  decouples  the  amplitude  of 
the  angular  profiles  from  their  shape,  through  the  use  of  Y  and  E*  respectively.  More  recent  work  by  Zhang  et  al28 
also  discusses  the  use  of  a  varying  energy  parameter,  but  in  the  context  of  expressions  for  energy  distributions  and 
total  sputter  yields.  In  general,  rather  than  using  the  MZ  expressions  for  a  priori  calculation,  we  treat  Y  and  E*  as 
free  fit-parameters  which  we  determine  from  (least- squares  fitting)  experimental  data.  Note  that  profile  shapes  are 
determined  by  the  ratio  EVE  and  for  high  ion  energy  ( EVE  «  1)  the  MZ  expression  reduces  to  the  diffuse  yield 
(y=Tcos(d)/7t). 
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F.  Weight  Loss  Analysis 


A  microgram  scale  was  used  to  weigh  the  samples  (before  and  after  each  sputtering  session).  With  this 
method  the  volumetric  sputter  yield,  Y  (in  units  of  mm3/C),  of  each  sample  is  found  using  Eqn.  (3),  where  Am 
represents  the  mass  loss  due  to  ion  beam  exposure,  t  is  the  exposure  time,  and  other  symbols  are  defined  as  above: 


Y  = 


Am 

P  J B,Avg  t 


(3) 


Measurement  uncertainties  for  total  sputter  yields  were  calculated  by  estimating  individual  contributions  and 
combining  them.  Uncertainties  are  as  follows:  ±  15%  on  beam  current,  ±  200  (Ltg  on  mass  measurements,  and  +1% 
on  test  time  lengths.  The  error  on  mass  tends  to  dominate  at  low  yields,  and  we  comment  that  the  quoted  mass 
uncertainty  is  an  order  of  magnitude  larger  than  that  provided  by  the  scales  themselves.  The  resultant  uncertainty 
varies  depending  on  measurement  conditions  but  typical  uncertainties  are  approximately  25%.  Note  that  this  value  is 
a  random  error  resulting  from  finite  accuracy  of  the  mass,  current,  and  time  measurements.  The  influences  of  other 
factors  affecting  the  accuracy,  such  as  humidity  variation,  are  not  taken  into  account  in  this  calculation.  See  a 
discussion  of  possible  factors  affecting  the  accuracy  of  the  weight  loss  measurements  in  Section  IV A. 


G.  Measurement  Procedure 

Targets  are  deliberately  pre-sputtered  to  better  represent  the  conditions  found  in  long-duration  EP  operating 
applications.  Pre- sputtering  for  HBC  and  HBR  targets  is  an  inherent  consequence  of  reusing  the  same  targets  used 
in  our  past  sputter  tests.  Prior  to  initial  testing,  the  (new)  HP  BN  target  was  pre-sputtered  to  clean  its  surface  and 
remove  any  potential  surface  layers.  Pre- sputtering  was  accomplished  by  a  750  eV  ion  beam  with  a  current  density 
of  ~1  mA/cm2  for  2  hours.  An  order-of-magnitude  estimate  for  the  typical  dose  of  incident  ions  on  a  target  prior  to 
testing  is  approximately  1020‘21  ions/cm2  (corresponding  to  10s- 100s  of  hours  and  eroded  thickness  of  approximately 
10-100  microns).  Target  contamination  effects  are  estimated  to  be  negligible,  since  for  typical  conditions  the  flux  of 
ions  incident  on  the  target  is  approximately  10  times  higher  than  the  flux  of  nitrogen  (the  major  contaminant)  to  the 
target7. 

Before  placement  into  the  vacuum  chamber,  targets  are  placed  on  the  scales  described  above.  Each  target  is 
weighed  several  times  to  ensure  stability,  and  the  average  value  of  the  measured  masses  is  taken.  Targets  are  then 
inserted  into  the  vacuum  chamber  and  the  sputter  measurement  is  conducted.  Test  durations  are  fixed  such  that  the 
mass  change  of  the  test  is  large  compared  to  the  precision  of  mass  measurements,  and  such  that  the  QCM  has  time 
(at  each  position)  to  sufficiently  stabilize  relative  to  thermal  and  background  noise.  Test  times  vary  from  several 
hours  for  the  higher  energy  tests  (250  eV)  to  as  long  as  14  hours  for  the  lower  energy  tests  (60  eV).  After  each  test, 
the  targets  are  again  weighed  and  the  difference  in  masses  is  used  for  the  weight-loss  measurement.  As  described  in 
Section  IIIC,  we  do  not  perform  weight  loss  measurements  for  the  HP  samples,  and  special  considerations  are 
employed  for  the  HBR  samples. 


IV.  Results  &  Discussion 

A.  Total  Sputter  Yields 

Validation  of  our  total  sputter  yield  measurements  was  performed  by  using  molybdenum  as  a  control.  While 
there  is  variation  in  data  from  different  research  groups,  molybdenum  sputter  yields  are  reasonably  well 
characterized36.  We  have  measured  the  total  sputter  yield  of  Mo  at  normal  incidence  at  energies  of  80,  150,  and  200 
eV.  The  total  yields  from  weight  loss  and  QCM  (found  from  the  Y  parameter  of  the  best-fit  MZ  profile)  were  self- 
consistent  and  in  agreement  with  the  Yamamura  and  Tawara  curve  fit36  to  within  experimental  uncertainty  (see 
Section  HIE).  These  agreements  are  taken  as  validation  of  the  measurement  methods  and  in  particular  of  the 
assumption  that  the  sticking  coefficient  of  the  QCM  may  be  taken  as  unity  for  condensables.  Figure  3  shows  total 
sputter  yields  measured  by  weight  loss  as  a  function  of  ion  energy  for  HBC  grade  and  HBR  grade  boron  nitride  for 
normal  incidence. 

It  should  be  mentioned  that  the  measurements  of  total  sputter  yield  values  were  performed  multiple  times  for  the 
majority  of  the  data  points.  The  reproducibility  of  the  QCM  results  was  good,  while  the  weight  loss  based  sputter 
yield  values  varied  significantly.  For  the  cases  when  several  measurements  were  done  at  the  same  conditions,  the 
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average  value  is  reported.  Random  measurement  errors,  discussed  in  section  IIIF,  cannot  explain  these  variations,  so 
the  error  bars  on  the  weight  loss  sputter  yield  plots  are  based  on  the  reproducibility  of  the  measurements,  rather  than 
on  the  values  described  in  section  IIIF.  The  mass  error  is  attributed  primarily  to  effects  of  humidity  and  moisture 
absorption.  The  work  is  ongoing  to  understand  moisture  pickup  effects  and  to  modify  the  measurement  procedure  in 
order  to  minimize  their  influence. 


E  (eV) 


Figure  3.  Total  sputter  yield  from  weight  loss  versus  ion  energy  for  normal  incidence.  100  eV  data  points 

coincide. 

Using  the  Y  parameter  of  the  best- fit  MZ  profile  (equivalent  to  integrating  the  best-fit  profile)  we  also  obtain 
total  sputter  yield  results  from  the  QCM  measurements.  Total  yields  found  in  this  way,  for  normal  incidence  ions, 
are  plotted  in  Fig.  4.  Recall  that  the  QCM  measures  only  condensable  components  that  stick  to  the  QCM,  which  for 
BN  corresponds  to  atomic  boron  and  BxNy  clusters  (including  BN),  but  not  nitrogen.  Therefore,  one  expects  total 
yields  from  QCM  measurement  to  be  below  the  weight  loss  values.  (Note  that  the  QCM  measurement  would 
actually  be  more  useful  in  certain  cases,  such  as  for  modeling  of  deposition  and  sticking  of  sputtered  particles  to 
spacecraft  surfaces.)  Following  these  ideas,  Fig.  5  plots  the  ratio  of  total  yield  from  weight  loss  to  the  corresponding 
yield  from  QCM  measurements.  We  find  that  the  ratio  of  the  weight  loss  to  QCM  yields  tends  to  be  >1.75  which 
corresponds  to  BN  sputter  products  being  predominantly  atoms.  If  all  sputtering  was  as  BxNy  clusters,  the  ratio  of 
yields  would  be  unity,  while  if  all  sputtering  was  as  atomic  boron  and  nitrogen  the  ratio  would  be  2.30  (= 
( Mb+Mn)IMb ,  where  MB  and  MN  are  atomic  masses  of  B  and  N  respectively).  Sputtering  as  atoms  is  consistent  with 
past  mass  spectrometry  results1  and  multi-component  sputtering  theory  suggests  sputtering  predominantly  as  atoms. 


E  (eV) 


Figure  4.  Total  sputter  yield  from  QCM  versus  ion  energy  for  normal  incidence. 
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Figure  5.  Total  sputter  yield  from  weight  loss  divided  by  total  sputter  yield  from  QCM  (see  text). 

We  have  studied  the  variation  of  total  sputter  yield  with  incidence  angle  at  a  fixed  energy  of  100  eV. 
Measurements  have  been  performed  with  both  weight  loss  and  QCM  and  results  are  shown  in  Fig.  6  and  Fig.  7 
respectively.  The  weight  loss  results  are  higher  than  QCM  results  for  the  reasons  described  in  connection  with 
Fig.  5.  For  the  QCM  data  the  trend  is  for  the  yields  to  somewhat  increase  over  the  range  of  angles  measured 
consistent  with  past  research37.  The  increase  with  angle  is  rather  weak  for  weight  loss  measurements.  This  is 
contrary  to  the  QCM  data,  to  our  past  measurements  and  to  other  published  data,  where  total  sputter  yield 
increases  with  the  incidence  angle.  The  reasons  of  this  discrepancy  are  not  clear  and  will  require  further 
investigation.  One  possible  explanation  is  that  due  to  the  graphite-like  structure  of  hot-pressed  BN,  the 
probability  of  xenon  ion  implantation,  i.e.  of  xenon  ion  becoming  stuck  between  the  molecular  planes  of  BN, 
increases  with  incidence  angle.  Xenon  ion  implantation  has  been  observed  in  the  past  during  sputter  yield 
measurements  of  ceramic  materials38.  If  the  number  of  xenon  atoms  stuck  in  the  target  increases  with  the 
incidence  angle,  then  while  the  actual  sputter  yield  also  increases,  the  superposition  of  the  two  trends  can  result 
in  almost  flat  measured  sputter  yield  vs.  incidence  angle  dependence.  Note  that  the  QCM  measurements  should 
not  be  affected  be  xenon  implantation,  which  is  in  agreement  with  the  data  presented  on  Fig. 7,  where  total 
sputter  yield  increases  with  the  increase  of  incidence  angle. 
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Figure  6.  Total  sputter  yield  from  weight  loss  versus  incidence  angle  for  100  eV  ion  energy. 
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Figure  7.  Total  sputter  yield  from  QCM  versus  incidence  angle  for  100  eV  ion  energy. 


Where  possible,  we  compare  our  measured  HBC  BN  (weight  loss)  total  sputter  yields  with  measurements  and 
modeled  values  from  other  research  groups.  Figure  8  shows  our  recent  measurements  along  with  weight  loss 

39  26  40  21 

measurements  by  Semenov  ,  Gamier  ,  Abashkin  ,  and  Yim  .  Note  that  not  all  of  the  authors  specify  the  grade  of 
BN  used.  We  note  that  the  values  we  are  currently  reporting  are  somewhat  higher  relative  to  past  measurements. 
Ongoing  work  is  examining  possible  systematic  errors  including  humidity  effects  mentioned  above. 
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Figure  8.  Total  BN  sputter  yields  as  compared  with  published  values. 

We  also  discuss  our  total  yield  measurements  in  comparison  to  past  measurements  from  our  laboratory9'11. 
We  find  that  the  currently  reported  total  yields  from  the  QCM  are  generally  within  error  of  the  previously 
reported  values.  However,  the  currently  reported  weight  loss  values  are  higher  by  a  factor  of  2-3  relative  to 
previously  reported  values.  The  reasons  for  this  discrepancy  are  unclear;  however,  we  do  comment  that  in  the 
past,  the  QCM  and  weight  loss  measurements  have  been  performed  in  different  chambers  and  not  concurrently. 
Further,  the  past  weight  loss  measurements  used  different  neutralization  schemes.  The  improvements  in 
experimental  configuration  and  consistency  of  data  suggest  that  that  the  current  results  should  be  considered  to 
be  more  accurate.  On  the  other  hand,  as  discussed  above,  the  values  appear  elevated  to  those  in  the  literature  and 
are  therefore  being  more  closely  examined. 
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B.  Differential  Sputter  Yields 


Using  the  QCM  measurement  technique  we  have  measured  differential  sputter  yield  profiles.  As  described  in 
Section  HIE,  the  MZ  expressions  are  used  to  fit  the  profile  shapes  and  the  parameter  EVE  describes  the  shape  of  the 
profile.  Figure  9  shows  the  variation  in  EVE  for  normal  incidence  as  a  function  of  ion  energy  for  the  three  grades  of 
BN.  Figure  10  shows  the  variation  in  EVE  at  100  eV  as  a  function  of  incidence  angle  for  the  three  grades  of  BN. 


Figure  9.  EVE  versus  ion  energy  for  normal  incidence. 


Figure  10.  EVE  versus  incidence  angle  for  100  eV  ion  energy. 

Examples  of  comparison  between  measured  (raw)  QCM  data  and  fitted  MZ  profiles  are  given  in  Fig.  11.  Both 
plots  are  for  xenon  ion  energies  of  100  eV  on  HBC  BN,  with  the  left  plot  being  at  normal  incidence  and  the  right 
plot  being  at  30  degrees  incidence.  The  plots  include  QCM  measured  points,  best- fit  MZ  profiles,  and  (for 
comparison)  diffuse  profiles  with  the  same  total  yield.  One  can  see  relatively  good  agreement  between  the 
measured  profiles  and  MZ  profiles.  The  normally  incident  profile  is  azimuthally  symmetric.  The  profile  for  30 
degrees  incidence  is  measured  in  the  forward/backward  plane  (0=0,180  degrees)  and  shows  a  forward  sputter  lobe 
(negative  alpha)  and  reduced  sputtering  in  the  backward  direction  (positive  alpha).  In  general,  the  MZ  expressions 
provide  reasonable  descriptions  of  the  measured  profiles,  but  tend  to  predict  a  slightly  broader  and  lower  amplitude 
forward-sputter  lobe. 
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Figure  11.  Example  of  QCM  data  with  best-fit  MZ  profiles  for  100  eV  ions  on  HBC  BN. 

Left:  Normal  incidence.  Right:  30  degrees  incidence. 

In  Fig.  12  we  plot  examples  of  the  best-fit  MZ  differential  sputter  yields  using  colored  hemispheres.  Colors 
(indicated  in  legend)  correspond  to  the  yield  in  the  given  direction.  For  normal  incidence  the  profile  is  azimuthally 
symmetric,  while  for  increasingly  non-normal  incidence  the  profiles  show  increasing  forward  sputter  lobes. 
Interestingly,  even  for  15  degrees  incidence,  the  profile  is  already  significantly  azimuthally  asymmetric.  Our  past 
work  has  measured  differential  sputter  yield  profiles  for  some  of  the  same  conditions10.  The  shapes  measured  here 
are  similar  to  those  of  the  past  work  with  EVE  ratios  consistent  to  approximately  +/-  0.2.  This  degree  of  consistency 
corresponds  to  qualitatively  similar  shapes  in  our  current  versus  past  measurements 


Figure  12.  Differential  sputter  yield  profiles  for  normal  incidence  at  different  beam  energies.  Left:  HBC- 

grade  BN.  Right:  HBR-grade  BN 
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Figure  13.  Differential  sputter  yield  profiles  for  100  eV  ion  energy  at  varying  incidence  angle.  Left:  HBC- 

grade  BN.  Right:  HBR-grade  BN. 

V.  Conclusions 

We  have  detailed  the  development  of  an  experimental  apparatus  targeted  at  low  energy  sputter  yield 
measurements.  We  have  reported  total  and  differential  sputter  yield  measurements  for  three  grades  of  BN:  HBC, 
HBR,  and  HP.  Measurements  have  been  performed  using  a  weight  loss  approach  (for  total  sputter  yield)  and  a 
QCM  measurement  approach  for  differential  sputter  yield  profiles  and  total  sputter  yield  of  condensable 
components.  Using  a  novel  four-grid  source  we  have  performed  the  first,  to  our  knowledge,  measurements  of  BN 
sputtering  in  the  100  eV  range  and  below.  Effects  of  humidity  (moisture  pickup)  and  xenon  ion  implantation  are 
possibly  affecting  the  accuracy  of  our  measurements  and  ongoing  work  seeks  to  mitigate  these  effects.  The 
experimental  error  on  our  current  measurements  is  rather  large  and  limits  our  ability  to  compare  with  published 
results.  We  do  observe  higher  yields  relative  to  published  data  and  the  reasons  for  this  are  under  ongoing 
investigation.  In  comparison  to  a  representative  refractory  metal  such  as  molybdenum,  we  find  comparable 
volumetric  yields  for  BN,  though  the  corresponding  BN  mass-  (or  atomic-)  based  yields  are  still  ~5x  lower  due  to 
the  lower  density  of  BN. 

Angular  dependence  of  sputtering  has  been  studied  using  the  QCM  measurement  apparatus  and  we  find  that 
the  MZ  profiles  provide  a  reasonable  description  of  the  measured  profiles.  The  shapes  (EVE  values)  for  the  three 
grades  of  BN  are  all  relatively  similar  and  show  azimuthally  symmetric  behavior  at  normal  incidence  and 
forward/backward  sputtering  features  at  oblique  incidence.  Comparison  of  total  yield  results  from  the  QCM  and 
weight  loss  measurement  systems  give  information  on  the  composition  of  sputtered  particles,  since  the  former 
method  captures  only  condensable  particles.  We  find  that  the  total  yields  from  weight  loss  are  roughly  twice  those 
of  the  QCM,  indicating  that  sputtering  is  largely  as  atoms  (with  boron  atoms  captured  by  the  QCM,  but  not  nitrogen 
atoms).  This  information  is  particularly  relevant  for  interpretation  of  future  laser  diagnostics  systems  based  on 
CRDS  which  will  measure  only  boron  contributions41. 

There  is  a  critical  need  in  the  EP  community  for  low  ion  energy  sputter  measurements  of  BN,  and  the  present 
contribution  is  a  step  forward  filling  this  gap.  The  results  presented  here  are  part  of  a  comprehensive  study  currently 
underway.  Upcoming  and  ongoing  work  aims  to  measure  the  BN  sputtering  over  a  broader  range  of  sputtering 
conditions  and  includes  study  of  variation  of  sputter  yields  with  BN  temperature. 
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Sputter  Erosion  Measurements 


•  Sputtering  is  key  mechanism  that  places  fundamental 
limitations  on  ion  and  plasma  thruster  life  times 

•  Sputtered  material  can  cause  space  craft  contamination 

->  Device  measurements  and  material  characterization 
measurements  needed  for  code  validation/development 

•  Available  Techniques:  Weight  loss,  deposition  monitor, 
surface  layer  activation,  profilometry,  LIF,  optical  emission, 

CRDS 
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BN  Erosion  in  SPTs 


•  Erosion  of  BN  insulator  leads  to  end  of  life  in  stationary 
plasma  thrusters  (SPTs). 

•  Low  energies  (~<1 00  eV  range)  critical. 
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Total  and  Differential  Sputter  Yields 


•  Total  yield: 

Y  (atoms/ion  ->  mm3/C) 

Y=  Y (E,  /?,  material  (grade),  7) 

where:  E  =  ion  energy,  j3=  incidence  angle,  T  =  temperature 


Differential  yield: 

V'  (atoms/ion/sr  ->  mm3/C/sr) 
y=Yx  y\a,<p) 

where:  a=  polar  ejection  angle,  </>=  azimuthal  ejection  angle 

2717112 

Y  =  }  y(a, 0)dQ.  =  f  f  y(a,<t>)sm(a)dad<f> 


0  0 


Diffuse: 

Or: 


y  =— cos(or) 

TC 

y  occosn(a) 
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Available  BN  Sputter  Data 


Limited  data  available,  especially  at  low  energy  (<~100  eV) 

•  Effect  of  grades? 

•  Temperature? 

•  Incidence  angle? 
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Sputter  Measurement  Apparatus 
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Sputter  Measurement  Apparatus 
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Sputter  Measurement  Apparatus 


•  Stainless  steel  vacuum  chamber  0.125  m3.  Baseline  pressure  below 
1*1  O'6  Torr 

•  CTI-8  cryogenic  pump  (1 500  liter/s  for  air) 

•  Hot  filament  ion  source  (four  grid  ion  optics  system) 

•  Neutralizers:  filament  -  ion  beam  neutralization,  and 
PBN  -  sample  surface  neutralization 

•  Sigma  Instruments  SQC-339  deposition  controller 
(reads  the  crystal  frequency  to  0.001  Hz) 

•  RC-cut  QCM  (Tangydine  corp.) 

•  Polyscience  9000  series  digital  temperature  controller 
that  allows  temperature  control  to  better  than  0.01  K 

•  Compumotor  motion  control  system 

•  Lab  View  based  data  acquisition  system 
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Four-Grid  Ion  Source  -  make  slide 


•  Our  old  (two-grid  system)  was  limited  by  -100  eV 

•  A  new  four-grid  ion  optics  enables  to  extract  collimated 
beams  at  the  low  ion  energies  (below  100  eV) 

•We  were  able  to  extract  well-collimated  beams 
•even  at  energies  as  low  as  30  eV 


4-Grid  Ion  Source  Beam  Current 
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QCM  Deposition  Monitor 
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Measurement  Procedure 


•  Mount  6  inch  square  samples  in  sample  holder 

•  Pump  down  to  high  vacuum  (10'6  Torr  range) 

•  Pre  sputter  to  remove  surface  layer  affected  by  polishing 

•  Vent  and  wait  for  the  sample  weight  to  equilibrate  (if 
necessary) 

•  Initial  weight  measurement 

•  Mount  6  inch  square  samples  in  sample  holder 

•  Pump  down  to  high  vacuum  (10'6  Torr  range) 

•  Run  ion  source  and  perform  QCM  differential  sputter  yield 
measurements 

•  Vent  and  wait  for  the  sample  weight  to  equilibrate  (if 
necessary) 

•  Weight  measurements 
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Materials  and  Conditions 


•  BN  is  material  of  primary  interest  for  Hall  thrusters  from 
the  point  of  view  of  erosion  and  contamination 

•  Boron  Nitride:  HBC  and  HBR  grades  (  GE’s  Advanced 
Ceramics),  HP  grade  (Saint-Gobain): 

-  HBR,  HP  -  calcium  borate  is  used  as  binder 

-  HBC  -  no  binder 

•  All  three  BN  grades  are  hot-pressed  graphite-like 
allotropes  of  BN.  Three  grades  have  similar  properties 
with  some  differences,  e.g.  HBR  and  HP  have  higher 
moisture  absorption  than  HBC 

•  Measurement  conditions:  Xe  ion  energies  of  60,  80,  100, 
200  and  250  eV,  incidence  angles  of  0, 15,  30,  and  45 
degrees 
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Sputter  Yield  Measurement  -  Weight  Loss 


Y  (mm3/C ): 


where 


Am  =  mass  change 

p  =  density  of  target  material 

JB  =  current  of  bombarding  particles 

t  =  exposure  time 
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Moisture  Effects:  HBC  and  HBR 


Hours  Hours 


Moisture  absorption  of  various  grades  of  hot-pressed  boron  nitride  (from  GE  Webpage). 
Humidity  in  the  lab  is  about  40%.  It  is  reasonable  to  expect  that  the  HBC  BN  samples  will  have 
essentially  no  moisture  pick  up,  but  the  HBR  samples  may  have  up  to  1%  pickup. 
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Moisture  Effects:  HP 


Time  (minutes) 


Mass  accumulation  as  a  function  of  time  from  six  different 
HP  BN  weighing  sessions.  Time  zero  corresponds  to  90 
minutes  after  removal  from  vacuum  chamber. 

Large  mass  buildup  and  variation  of  individual  samples  relative  to  the  total  before/after 
mass  change  of  a  typical  test  (-2-10  mg)  precludes  weight  loss  measurements  of  HP 
BN  unless  proper  measurement  procedure  is  found 
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Sputter  Yield  Measurement  -  QCM 


Y  (mm3/C/sr ): 


where: 


R(a<p,)  =  measured  mass  accumulation  rate 
p  =  density  of  target  material 
JB  =  ion  current  incident  on  the  target 
r  =  distance  from  target  center  to  QCM 
As  =  QCM  sensor  area  (0.535  cm2) 
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Profile  Fitting  -  Modified  Zhang  Eqns. 


y  =  differential  sputter  yield 
Y  =  total  sputter  yield  -  FIT  PARAMETER 
E*=  characteristic  energy  -  FIT  PARAMETER 


p=  incidence  angle 
a=  polar  ejection  angle 
</>=  azimuthal  ejection  angle 
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Y  (atoms/ion) 


Sputter  Measurement  Validation 


Total  sputter  yield  for  normally  incident  Xe+  on  Mo 
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Y  (mmA3/C) 


Total  Sputter  Yield  versus  Energy  (Weight  Loss) 


E  (eV) 


Normal  Incidence 


Distribution  A:  Approved  for  public  release;  distribution  unlimited 


0/evxuxu)^ 


Total  Sputter  Yield  versus  Energy  (QCM) 


E  (eV) 


Normal  Incidence 
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Y(WL)  /  Y  (QCM) 


Yi©ld  weight  Loss^  Yi6ld 
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{Mb+Mn)/Mb  =  2.30 
Mb  -  atomic  mass  of  B 
Mn- atomic  mass  of  N 
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0/£vUIUI)  X 


Total  Sputter  Yield  versus  Angle  (Weight  Loss) 


P  (degrees) 


E=  100  eV 
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Y  (mmA3/C) 


Total  Sputter  Yield  versus  Angle  (QCM) 


P  (degrees) 


E=  100  eV 
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0/£vUIUl)  A 


Comparison  with  Past  Measurements 


E(eV) 


Past  measurements  limited  to  minimum  of  ~100  eV. 
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Differential  Sputter  Yield  -  Xe+  on  HBC  BN 


E=  100  eV,  /?=  0° 
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Y  (mm3/C/sr) 


Differential  Sputter  Yield  -  Xe+  on  HBC  BN 


E  =  100  eV,  p=  30°  (Normal) 


Forward  Sputter  ^ 


Backward  Sputter 
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Xe+  on  HBC  BN  -  Differential  Yield  Profiles 


100  eV 
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EVE 


Differential  Yield  Profiles  -  EVE  Shape  Parameter 


E(eV) 


p  (degrees) 
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Conclusions 


•Use  of  four-grid  ion  source  with  weight-loss  and  QCM  allows  sensitive 
measurement  of  total  and  differential  sputter  yields  for  different  grades  of 
BN  at  ion  energies. 

•Total  sputter  yields  for  all  grades  are  similar. 

•Comparison  of  weight  loss  to  QCM  yields  indicates  sputtering  primarily  as 
atoms. 

•Total  sputter  yields  in  1 00-250  eV  range  are  somewhat  higher  relative  to 
those  in  the  literature. 

•Measurements  in  the  <100  eV  range  are  critical  to  thruster  erosion  and 
reported  measurements  are  amongst  the  first. 

•Differential  sputter  yield  profiles  for  different  grades  are  relatively  similar 
and  show  azimuthally  symmetric  behavior  at  normal  incidence  and 
forward/backward  sputtering  features  at  oblique  incidence. 
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Ongoing  and  Future  Work 


•  Expand  test-matrix  of  energies,  angles 

•  Validation  studies  (moisture,  neutralization 
effects,  Xe  ion  implantation) 

•  Improved  precision  in  mass  measurement 

•  Studies  of  dependence  of  sputter  yields  on  target 
temperature 
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Extra  Slides 
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Chamber  and  Pressures: 

•1500  liter/s  CTI-8  cryogenic  pump 

•Chamber  base  pressure  -  5><1 0-7  Torr  (after  8-hour  bake  out) 
•Working  pressure  -  0.6-1  xl  0"4  Torr 

Ion  Source  and  Beam: 

•  Typical  Xe  flowrate  -  0.3-0.5  seem 

•  Four-grid  ion  optic  system 

•  Thoriated  tungsten  hot-filament  cathode 

•  Discharge  voltage  ( VD)  -  30-40  V 

•  Average  current  densities  ~0.1  mA/cm2 

Sample  Preparation: 

•Sputter-cleaned  for  3-6  hours  with  a  500-750  eV  (~0.2  mA/cm2) 
beam 

PBN: 

•Emission  current  -  8-20  mA 
•Xe  mass  flow  rate  =  0.5  seem 
•Voltage  bias  =  -24  V 
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Azimuthally  Symmetric  Profiles:  y(a) 


•  Simplest  case:  diffuse  (cosine)  profile.  (Rare,  but  convenient) 

•  Over-cosine  profiles:  high  energy  ions  (all  incidence  angles) 

•  Under-cosine  profiles:  low  energy  ions  (normal  incidence) 


-0.6  -0.4  -0.2  0  0.2  0.4  0.6 

y  (atoms/ion/steradian) 
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Azimuthally  Symmetric  Profiles:  y(a) 


350  eV  Xe  ions  normally  incident  on  Mo 
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Non-Azimuthally  Symmetric  Profiles:  y(a,</>) 


Anisotropic  targets 

Low  energy  sputtering  (oblique  incidence  ions): 
•  Isotropic  sputtering  condition:  (M/E//Wf£/)1/2«  1 


350  eV  Xe  ions  30  degree  incident  on  Mo 


Back 
Scatter 
Direction 
^=180a 

Distribution  A:  Approved  for  public  release;  distribution  unlimited 


Coliseum  Expressions  (for  reference) 


Constant: 


Jc 


(E,0)  =  c 


Kannenberg: 


yK 


(E,  e) = (c, + ctiE)[c,o + CQ02 + C;+ ) 


Pencil: 


Roussel:  JR (e,  J)  =  Ic~A  +  C^X1  - 0.720 + 1 1  ,7202  - 3. 1 303  - 2.5704 ) 


■Yamamura: 


yr 


{E,0)  =  { 


/□  -l(cos  1 0- 1 

cos  7  tfe  v 


1  +  0.35J7  i  (f) 


1- 

( E  \ 

^ th 

0.5 

E  j 

-|2.8 


/=/, 


i-C 


1  +  2.5 

<r 


L  =  f  cos  0t 


opt 


^  _  j  _  I  E th,A 


J 
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Total  Yield  Expressions  (Normal  Incidence)  -  Y&T  2/5 


•  Example  of  comparison  of  Y&T  expressions  to  CSU  and  other  data: 


Ion  Energy,  eV 


Figure  I  I.  Total  viclcls  for  Xe  on  Ta. 


From:  Differential  Sputtering  Yields  of  Refractory  Metals  by  Xenon,  Krypton,  and  Argon  Ion  Bombardment 
at  Normal  and  Oblique  Incidences,  IEPC-2005-293,  Zoerb  et  al,  Colorado  State  University 
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Total  Yield  Expressions  (Varied  Incidence) 


•  Normal  incidence  total  yields  Y(0=0)  can  be  scaled  to  give  total  yields  Y(0)  for 
varying  angle  of  incidence  6. 

•  Use  approach  of  Yamamura  (as  shown  on  the  Coliseum  slide).  Simple  scaling 
expression  is  based  on  tabulated  values  of  T  and  “<9op/\ 

See:  Y.  Yamamura,  Y.  Itikawa,  N.  Itoh,  “Angular  Dependence  of  Sputtering 
Yields  of  Monatomic  Solids”  IPPJ-AM-26  (1983) 

•  In  this  case  the  database  of  material  constants  is  relatively  limited.  For  other 
ion/target  materials,  the  (2)  needed  parameters  can  be  easily  found  by  fitting 
CSU  (or  other)  data  to  these  functional  forms. 


•  Expression  and  example: 


y(e) 

Y(0) 


f-  cos  (©opt) 


1 

cos(g) 


Example  is  Ar  on  Cu  at  500  eV  (table  3  of  paper) 

f  :=  3.35  0opt  :=  49.6 


Y  peaks  at  #=#opt  (=49.6° 
for  Ar+  on  Cu  at  500  eV). 
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Differential  Sputter  Yield  Expressions  -  3/7 


•  Base  expressions  are  from: 

ANISOTROPIC  ANGULAR  DISTRIBUTION  OF 
SPUTTERED  ATOMS 

ZHU  LIN  ZHANG11*  and  LAI  ZHANGb 

a  Department  of  Mathematics  and  Physics  and  b Department  of  Computer  Science  and  Technology, 
Anhui  University  of  Science  and  Technology,  Huainan,  China 

•  Differential  sputter  yield,  y(a,</>): 


,  ^  n^a(M2/M')sAE)  ,  , 

y  ( a,  (p)  -  0.042  — - - - cos  ( a) 


7lU  „ 


i-^^fcost^rfaJ  +  l^sinfAfsin 


where 


.  v  3sin(a)2-l  cos(a)  (3sin(ar)  +lj  fl  +  sin (a) 
y{a)  = -  a —  + - - - i - -In 


sin(a)^ 


2  sin  (a) 


V 


l-sin(a) 


A 


J 
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Differential  Sputter  Yield  Expressions  -  4/7 


Problem  1 :  Total  sputter  yields 

•  The  total  sputter  yields  in  the  Zhang  expressions  (based  on  earlier  work  by 
Yamamura  and  others)  do  not  have  the  flexibility  to  use  free  (fit)  parameters 
and  do  not  agree  well  with  data.  The  total  yields  are  not  explicitly  included,  but 
can  be  determined  as: 


l7l  7lH 


Y(E,0)=  f  f  y(a,</>) sin (a)dad(/>  =  0.042 


a(M2IM1)s,{E) 


0  0 


u. 


1- 


Jth 


COS 


(A) 


This  can  be  compared  with  the  more  complex  form  of  Y&T. 

Example  of  disagreement  of  “simple”  Zhang  total  sputter  yield  with  “better” 
Y&T  total  sputter  yield,  computed  for  Xe+  on  Mo: 


Note:  I  have  used  the  Y&T  value  for  Eth  to 
plot  the  Zhang  total  yield.  Zhang  suggests 
fitting  his  curve  to  data  to  obtain  values  of 
Eth,  but  (at  least  in  this  case)  such  an 
approach  also  does  not  allow  good 
agreement  with  the  Y&T  yields  (which  do 
agree  well  with  measurement). 


10  100  1 -103  1 -104 


Ion  Energy  (eV) 
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Differential  Sputter  Yield  Expressions  -  6/7 


Problem  2  -  Threshold  Energy 

•  Comparison  of  experimental  differential  sputter  yield  profiles  with  those 
calculated  from  Zhang  expressions  using  “a  priori’  Eth  values  (from  Y&T)  does 
not  give  good  agreement. 

•  Example:  Xe+  on  Mo,  45°  incidence,  750  eV 


Experimental  (Polar)  Surface  Zhang  (Polar)  Surface  using  “a  priori'  Eth= 46.83  eV 


Notes: 

1 .  Experimental  surface  from  CSU  Quartz  Crystal  Microbalance  Data  (see  AIAA  4336-2006). 

2.  Surfaces  are  plotted  with  same  scales  and  correspond  to  same  total  yield  (total  yield  is  not «  volume). 

3.  Zhang  surface  is  too  “diffuse”  owing  to  large  E/Eth  ratio  (see  next  slide). 
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Differential  Sputter  Yield  Expressions  -  7/7 


Solution  of  Problem  2  -  Fit  experimental  profiles  to  find  E*  values 

•  Improved  agreement  between  measured  angular  profiles  and  the  Zhang 

expressions  can  be  found  by  using  “fitted”  values  of  E*  in  place  of  Eth.  (We  do 
a  least  squares  fit  to  find  value  of  E*  giving  best  agreement  between 
experimental  profile  and  y(Y,E,E*,fi,a,</)).  For  a  given  ion  target  combination 
we  find  that  E*  varies  with  ion  energy  (likely  due  to  surface  texturing  effects, 
possible  variation  of  binding  energy  with  incident  ion  energy  etc.) 

Experimental  (Polar)  Surface  Zhang  (Polar)  Surface  using  fitted  E*  =  390  eV 


Notes: 

1 .  Zhang  surface  with  E*  better  captures  forward  scatter  (non-diffuse)  behavior  (though  is  still  slightly  too 
diffuse,  see  AIAA  4336-2006).  For  cases  studied,  agreement  between  profiles  improves  to  <std><~1 5%. 

2.  Both  surfaces  are  plotted  with  same  scales  and  correspond  to  same  total  yield. 
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